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ABSTRACT 

We present new IRAM 30m spectroscopic observations of the ~ 88 GHz band, including emission 
from the GGH (N = 1 —>• 0) multiplet, HGN (J = 1 —0), HCO+ (J = 1 —>■ 0), and HNG (J = 1 —>■ 0), 
for a sample of 58 local luminous and ultraluminous infrared galaxies from the Great Observato¬ 
ries All-sky LIRG Survey (GOALS). By combining our new IRAM data with literature data and 
Spitzer/IRS spectroscopy, we study the correspondence between these putative tracers of dense gas 
and the relative contribution of active galactic nuclei (AGN) and star formation to the mid-infrared 
luminosity of each system. We find the HGN (1-0) emission to be enhanced in AGN-dominated sys¬ 
tems (( ) = 1.84), compared to composite and starburst-dominated systems 

(( Lypj^ (i-o)/^HCo+ ( 1 - 0 ) ^ “ 1.14, and 0.88, respectively). However, some composite and starburst 
systems have L(jgj,j (i-o)/^hco+ (i-o) comparable to those of AGN, indicating that enhanced 

HGN emission is not uniquely associated with energetically dominant AGN. After removing AGN- 
dominated systems from the sample, we find a linear relationship (within the uncertainties) between 
logj^Q(L(j(^j^ ( 1 - 0 )) logiQ(LiR), consistent with most previous findings. L(jqj,^ typically in¬ 

terpreted as the dense gas depletion time, appears to have no systematic trend with Ljr for our sample 
of luminous and ultraluminous infrared galaxies, and has significant scatter. The galaxy-integrated 
HGN (1-0) and HCO"'' (1-0) emission do not appear to have a simple interpretation, in terms of the 
AGN dominance or the star formation rate, and are likely determined by multiple processes, including 
density and radiative effects. 

Subject headings: galaxies: ISM — galaxies: starburst — galaxies: active 
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Molecular gas is observationally linked to ongoing star 
formation through observed correlations between the star 
formation rate surface density and the H 2 surface den¬ 
sity as inferred from GO observations (e.g., Bigiel et al. 
2008; Leroy et al. 2012). GO (1-0) has a relatively low 
critical density {ncrit ~ 2 x 10^ cm“^) and so traces the 
bulk of the molecular gas. Molecular transitions such 
as HCN (1-0) and HCO'*' (1-0) have critical densities 
Hcrit ~ 3 X 10® cm“® and 2 x 10® cm“®, respectively, 
at 30 K, and so they are associated with higher density 
molecular hydrogen. Early studies found a linear corre¬ 
lation between HCN (1-0) and the infrared luminosity— 
Lir[8 — 1000 /rm]—of galaxies (Solomon et al. 1992; Gao 
& Solomon 2004b). This relation, which is tighter than 
that for GO (1-0) with Lir, was interpreted as evidence 
that HGN (1-0) traces the dense gas directly associ¬ 
ated with star formation. Revisiting the relationship 
between SFR and this molecular tracer, Garcia-Burillo 
et al. (2012) find Lrir to be a super-linear function of 
HGN (1-0), suggesting that there are other physical fac¬ 
tors that are important for the HCN emission. 

Calling into question the use of HCN (1-0) as a tracer 
of dense gas, several studies of systems hosting active 
galactic nuclei (AGN) have found integrated (Gracia- 
Garpio et al. 2006; Krips et al. 2008) and spatially- 
resolved enhancements (Kohno et al. 2003; Imanishi et al. 
2006, 2007, 2009; Davies et al. 2012) of HGN (1-0) emis¬ 
sion compared to what is observed in starburst galaxies. 
Similar results have been found for HGN (4-3) (Iman- 
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ishi & Nakanishi 2013, 2014). These results suggest that 
HCN emission is enhanced in the presence of AGN, po¬ 
tentially invalidating the use of HCN as a tracer of the 
dense molecular gas associated with ongoing star forma¬ 
tion, particularly in systems with AGN. However it is 
notable that some systems with a known AGN do not 
show this enhanced ratio (ARP 299, Imanishi & Nakan¬ 
ishi 2006; 1 Zw 1, Evans et al. 2006), suggesting the 
HCN (1-0) emission enhancement observed in some AGN 
hosts may have a more nuanced interpretation. 

Modeling by various authors suggests the HCN and 
HCO+ emission are affected by density, radiative, and 
abundance/ionization effects which potentially compli¬ 
cates interpretation of the line ratios in both photon 
dominated and X-ray dominated regions (PDR and 
XDR, respectively; e.g., Aalto et al. 1994, 1995; Huet- 
temeister et al. 1995; Lepp & Dalgarno 1996; Meijerink 
et al. 2007). The possible influence of the XDR on molec¬ 
ular abundances has been used to argue that elevated 
HCN (1-0) is a signpost of an AGN, however Juneau 
et al. (2009) and Costagliola et al. (2011) argue the exci¬ 
tation of HCN and HCO’*' is not solely driven by abun¬ 
dance and so this ratio may not solely trace XDRs. The 
relative influence of the other effects (infrared pumping, 
source compactness) has not been fully established. 

In order to investigate the relation between the pres¬ 
ence and strength of an AGN, and the HCN and HCO+ 
emission in a large sample of gas-rich galaxies, we have 
observed 58 luminous and ultraluminous infrared galax¬ 
ies ((U)LIRGs; Lir > 10^^ L 0 ), selected from the Great 
Observatories All-sky LIRG Survey (GOALS; Armus 
et al. 2009), with the Institut de Radioastronomie Mil- 
limetrique (IRAM) 30m Eight Mixer Receiver (EMIR; 
Carter et al. 2012). These new observations (Section 2) 
of the molecular rotational transitions, HCN (1-0) and 
HCO+ (1-0), are used in combination with calibrated 
AGN strengths determined from Spitzer/IRS (Houck 
et al. 2004) spectroscopy of the GOALS sample (Stier- 
walt et al. 2013; Inami et al. 2013) to assess the cor¬ 
relation of AGN strength with the global HCN (1-0) 
and HCO^ (1-0) emission (Section 3). We then dis¬ 
cuss the possible explanations for enhanced HCN (1-0) 
(Section 4). Finally, we explore the utility of HCN (1- 
0) and HCO+ (1-0) as tracers of the mass of dense gas 
associated with star formation (Section 5). 

The power of this study comes from the increased sam¬ 
ple size of (U)LIRGs with measurements of these lines 
and, particularly, from the large number of sources with 
measured mid-infrared diagnostic of the relative contri¬ 
butions of AGN and star formation to the infrared lumi¬ 
nosity of each system (a factor of 4 increase over Costagli¬ 
ola et al. 2011). The use of the mid-infrared diagnos¬ 
tics facilitates a good estimation of the importance of 
an AGN to the mid-infrared emission (i.e., as opposed 
to simply using rudimentary optical “AGN” and “star- 
burst” diagnostics to classify systems). This enables a 
direct investigation of the global 

line ratio (hereafter, HCN/HCO'’") as a function of the 
contribution of the AGN to the bolometric luminosity. 

Throughout the paper we adopt a WMAP-5 cosmol¬ 
ogy (Hinshaw et al. 2009, Hq = 70 km s“^ Mpc“^, 
f^vacuum = 0.72, = 0.28), with velocities cor¬ 

rected for the 3-attractor model of Mould et al. (2000). 


2. SAMPLE, OBSERVATIONS, AND DATA 

2.1. GOALS 

As noted, the data presented here were obtained as 
part of a millimeter survey of (U)LIRGs selected from 
GOALS. The GOALS sample as a whole consists of all 
the (U)LIRGs from the Revised Bright Galaxy Sample 
(RBGS; i.e., 60 fim flux density greater than 5.24 Jy 
and Lir > 10^^ Lq; Sanders et al. 2003). GOALS^® is 
a multi-wavelength survey aimed at understanding the 
physical conditions and activity in the most luminous 
galaxies in the local Universe. The dataset includes spec¬ 
troscopic and imaging observations in the infrared from 
Spitzer (Petrie et al. 2011; Inami et al. 2013; Stierwalt 
et al. 2013, J. M. Mazzarella et al. in prep) and Herschel 
(Diaz-Santos et al. 2013, 2014; Lu et al. 2014), GALEX 
and HST UV, optical, and near-infrared imaging (How¬ 
ell et al. 2010; Haan et al. 2011; Kim et al. 2013, A. 
S. Evans et al. in prep), Chandra X-ray observations 
(Iwasawa et al. 2011), and a suite of ground-based radio 
and sub-millimeter observations (e.g., Leroy et al. 2011; 
Barcos-Muhoz et al. 2015, Herrero-Illana et al. in prep). 
These observations collectively trace the obscured and 
unobscured activity and constrain the structural proper¬ 
ties and merger stages of these systems. 

The objects in this study were selected from two por¬ 
tions of GOALS; a high-luminosity sample (Lir > 10^^'^ 
L 0 ) to capture the most extreme star forming systems, 
and a sample of LIRGs with Lir < 10^^’^ L 0 which were 
selected to be isolated or non-interacting, on the basis 
of their stellar morphology (Stierwalt et al. 2013). The 
combination of these subsets of GOALS ensures this sam¬ 
ple spans the range of Lir and merger stage within the 
GOALS sample as a whole. 

2.2. Observations and Data Reduction 

The IRAM 30m Telescope was used with the EMIR 
receiver to observe an 8 GHz instantaneous bandwidth, 
tuned to simultaneously capture HCN (J = 1 —>• 0), 
HCO+ (J = 1 —>■ 0), the CCH (N = 1 —>• 0) multiplet, and 
HNC (J = 1 —0), thus reducing systematic uncertain¬ 
ties when comparing the fluxes of these lines. The rest 
frequencies^® for these transitions are: r'rest(HCN (1-0)) 
= 88.631 GHz, i^^e3t(HC0+ (1-0)) = 89.189 GHz, and 
i^r-est(HNC (1-0)) = 90.663 GHz. For CCH, we adopt a 
rest frequency which is the simple arithmetic mean of the 
individual multiplet frequencies, (CCH) = 87.370 
GHz. The beam FWHM for these measurements is 
^ 28”, corresponding to a linear size of 16.8 kpc at the 
mean redshift of our sample - thus we obtain galaxy- 
integrated measurements. Observations were done in 
wobbler switching mode (~ 1' throw, 0.8 s switching) to 
improve the baseline calibration. In Table 1 we list the 
pointing coordinates, assumed redshift, observing year & 
month, integration times, system temperatures, and the 
backend used, for each source. 

Some of our observations use the FTS backend, which 
consists of 24 individual fourier transform spectrometer 
units. In some scans, the FTS units exhibited gain vari¬ 
ations, resulting in offsets of the baselines for individual 

http;//goals.ipac.caltech.edu 

obtained from the Splatalogue database: 
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units. For those scans, linear fits were made to the base¬ 
line of each unit (masking out the expected locations of 
spectral lines), and the units were corrected to a common 
baseline. Scans were combined and exported using the 
GILDAS/CLASS software^^. Further analysis, includ¬ 
ing a linear baseline subtraction, smoothing, and flux 
measurement, was done using the pyspeckit software^® 
(Ginsburg & Mirocha 2011). The typical velocity reso¬ 
lution of these smoothed data is 60 MHz, corresponding 
to roughly 215 km s“^. In Figure 1 we show the reduced 
spectra for all sources, with the expected positions of the 
CCH, HGN (1-0), HGO+ (1-0), and HNC (1-0) lines 
marked. Lines were identified visually based on cata¬ 
loged optical redshifts and total fluxes were determined 
by integrating under the line. IRAS F16164-0746 had 
spectra which were clearly shifted from the location ex¬ 
pected based on the cataloged redshift. We measured a 
new redshift of 0.0229 for IRAS F16164-0746 compared 
to the value from NED, 0.0272. 

We determined uncertainties in the fluxes by calcu¬ 
lating the RMS, ach, per channel of width dv, in the 
line-free regions of the spectrum and multiplying by the 
square root of the number of channels, N^h covered by 
each line. We considered a line a detection if the mea¬ 
sured flux exceeded 3achdv^Nch- For non-detections we 
calculated a Scr upper limit by assuming a square line 
profile with a width half that of the detected HCN (1- 
0) or HCO+ (1-0) line, or 200 km s“^ if no lines were 
detected for that source. 

For IRAM 30m observations, temperatures are re¬ 
ported on the scale; these values were converted to 
a main beam temperature, T^b =T^/(Be///Fe//) us¬ 
ing the main beam efficiency, ^eff = 0.81, and a for¬ 
ward efhciency, Fg// = 0.95, both measured by IRAM 
staff on 26 August 2013^®. We use a Jy/K conversion of 
3.906(Fe///Ae//) = 6.185, where Ag// = 0.6, to convert 
the reported TJj values to flux densities. Line luminosi¬ 
ties were computed according to Solomon et al. (1992, 
their Equation 3), in units of K km s“^ pc^. 

As noted, the wide EMIR bandwidth enables simul¬ 
taneous measurements of HCN (1-0), HCO®' (1~0), 
HNC (1-0), and CCH. We had a detection rate of 78%, 
76%, 35%, and 37%, for HCN (1-0), HCO+ (1-0), 
HNC (1-0), and CCH, respectively. For detected lines 
the signal-to-noise weighted mean HNC (l-0)/HCN (1- 
0) and CCH/HCN (1-0) ratios are 0.5±0.3 and 0.8±0.3, 
respectively. One source, NGC 6285 has detected CCH 
emission and no detection of HCN (1-0) or HCO+ (1-0); 
the CCH detection is ~ 3.8 ct and, if confirmed, would be 
an interesting and rare example of a source with CCH 
brighter than HCN or HCO+. Such elevated CCH emis¬ 
sion may indicate an overabundance of CCH or a lack of 
dense molecular gas (Martin et al. 2014). The HCN (1- 
0)/HCO+ (1-0) ratio is discussed in more detail in Sec¬ 
tion 3. Integrated fluxes (or 3tT upper limits) are pro¬ 
vided for all lines, in Table 2, but we limit the analysis 
here to the HCN (1-0) and HCO®' (1-0) lines. Upper 
and lower limits in Figures use these 3(7 limits. CO (1- 
0) observations were also obtained as part of this pro¬ 
gram; these will be presented in R. Herrero-Illana et al. 

http://www.iram.fr/IRAMFR/GILDAS 

http://pyspeckit.bitbucket.org 

http://www.iram.es/IRAMES/mainWiki/IramSOmEfficiencies 


{in prep). 

In order to explore the influence of AGN on the 
HCN (1-0) and HCO®" (1-0) emission, we use the equiv¬ 
alent width (EQW) of the 6.2 p,m polycyclic aromatic 
hydrocarbon (PAH) measured from Spitzer/IRS low- 
resolution observations (Stierwalt et al. 2013). The EQW 
of the 6.2 pm PAH feature compares the strength of 
the PAH emission associated with star formation with 
the strength of the mid-infrared continuum (e.g., Genzel 
et al. 1998; Armus et al. 2004, 2007). The AGN is en¬ 
ergetically more important in systems with lower values 
of the PAH EQW. Points in Figures 2-7 are color-coded 
by their 6.2 pm PAH EQW, to show the relationship 
between AGN and starburst dominated systems. 

In addition to our new measurements, we incorporate 
observations from Gracia-Carpio et al. (2006), Costagli- 
ola et al. (2011), and Garcia-Burillo et al. (2012) in our 
analysis. The combination of these three samples in¬ 
cludes 63 individual galaxies with measured 6.2 pm PAH 
EQWs and detections in at least one of HCN (1-0) or 
HCO+ (1-0), comprising roughly 25% of the objects in 
the GOALS sample. 

We use Lir measurements from IRAS observations 
(Armus et al. 2009). For galaxies in pairs, we assign 
a portion of Ljr to each component based on their 70 
or 24 pm flux ratio, as described by Diaz-Santos et al. 
(2013). Measurements of the [C ii] 158 pm line are taken 
from Diaz-Santos et al. (2013). 

2 . 2 . 1 . Comparison with Previous Measurements 

Several of the sources observed by us have pre-existing 
published fluxes in the literature (e.g., Gao h Solomon 
2004a; Gracia-Carpio et al. 2006; Garcia-Burillo et al. 
2012). We compared our fluxes with those previous ef¬ 
forts and found good agreement for some sources (e.g., 
ARP 220, VV 114, NGC 0034, IRAS 15107-k0724, IRAS 
23365-1-3604), but we measure fluxes for some individual 
sources (NGC 2623, NGC 6701, UGC 5101, and NGC 
7591) which are factors of 2 — 3 higher than those previ¬ 
ously published. There appears to be no trend of these 
high fluxes with observing date or observing setup. De¬ 
spite some differences in fluxes for some sources, we find 
the HCN/HCO+ ratios for those sources are in good 
agreement with previously published ratios. 

2 . 3 . Distance and Aperture Effects 
2 . 3 . 1 . Source Distance Effects 

The 28" IRAM 30m beam FWHM covers linear scales 
of 6.1 — 46 kpc for these sources. Thus, these single-dish 
observations average together emission from many giant 
molecular clouds (GMCs) within each system, likely with 
varying physical conditions and environments. This fac¬ 
tor of ~ 8 in distance leads to possible concerns about the 
effects of beam size on the results, particularly if there 
are systematic variations in the CMC properties or en¬ 
vironments as a function of distance from the nuclei. 

Previous interferometric observations with ~ 5" reso¬ 
lution have found the HCN (1-0) and HCO®" (1-0) emis¬ 
sion to be unresolved in ULIRGs (Imanishi et al. 2007, 
2009). For the sources from their study that overlap 
with our sample, we find general agreement between the 
fluxes, suggesting the amount of extended flux is mini¬ 
mal and the HCN and HCO®" emission is confined to a 
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Figure 1. IRAM 30m EMIR spectra of the systems observed as part of this program, sorted by Right Ascension. In all panels we plot the 
observed as a function of the observed frequency. The typical velocity resolution of these smoothed data is 60 MHz, corresponding to 
roughly 215 km s“^. Dashed vertical lines mark (from left to right) the expected locations of CCH, HCN (1-0), HCO"*" (1-0), and HNC (1- 
0), based on optical redshifts from NED, except for IRAS F16164-0746, where we have used our new measured redshift. Information on 
the observing parameters are provided in Table 1. Integrated fluxes (or upper limits) for each line are given in Table 2. This figure is 
continued at the end of the manuscript, in the Figures currently labeled 8-10. 
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Figure 2. Ratio of the (1-0) ^HCO'^ (1 0) 

versus luminosity distance. The lack of a significant correlation 
suggests our results are not being affected by the distances to each 
system. Points are color-coded by the 6.2 fim PAH EQW. 

region smaller than the 28" beam of these IRAM 30m 
observations. 

In Figure 2 we compare the HCN/HGO’*" ratio with the 
luminosity distance; there is no obvious trend with dis¬ 
tance, suggesting that the projected beam size is not the 
dominant factor in determining the HCN/HCO'*' ratio in 
our sample. 

2 . 3 . 2 . Spitzer vs IRAM 30m Telescope Aperture Comparison 

The 6.2 pLm PAH EQW used to assess the relative dom¬ 
inance of the AGN in the mid-infrared was measured 
from a 3.6" aperture (Stierwalt et al. 2013). Thus, these 
measurements are upper limits to the large-scale influ¬ 
ence of the AGN, and we expect any systematic effect 
of this aperture difference between the mid-infrared and 
millimeter observations would serve to overestimate the 
importance of the AGN. 

In other words, we would expect any signature of AGN- 
dominated gas to be diluted with increasing source dis¬ 
tance. As will we shown (Section 3), the excitation 
of these molecular transitions does not appear to be 
solely a function of the AGN strength, thus we con¬ 
clude our results are not being significantly biased by 
the mis-match in aperture between the millimeter and 
mid-infrared datasets. 

3. ENHANCED GLOBAL HCN (1-0) EMISSION DOES NOT 
UNIQUELY TRACE AGN ACTIVITY 

Some of the initial claims that HCN (1-0) is 
enhanced in AGN were supported by plotting the 
^HCN (i-o)/^co ( 1 - 0 ) ^HCO+ (i-o)Tco ( 1 - 0 ) ^atios as 
a function of Lir. In Figure 3 we show the HCN (1- 
0)/HCO+ (1-0) ratio as a function of Lir^°, to compare 

Here we do not compare with CO (1-0) as this would po¬ 
tentially include molecular gas which is not physically associated 
with the regions emitting in HCN (1-0) and HCO+ (1-0). For this 
study, we rely on the HCN/HCO"*" ratio alone, to avoid concerns 
with mis-matched apertures between the CO (1—0) measurements 



Figure 3. Ratio of (1-0) (1 0) 

sus the total infrared luminosity. Points are color-coded by the 
6.2 pm PAH EQW. 

we previous work. As in previous studies (e.g., Gracia- 
Carpio et al. 2006), we find a relative dearth of sources 
at high Lir with low HCN/HCO+. 

However, our larger sample also contains a number of 
lower Lir systems with HCN/HCO’*' ratios as high as 
those of ULIRGs (> 1), which suggests that increasing 
Lir (or an associated increase in AGN contribution) may 
not be the sole driver of an increase in HCN/HCO"*' ra¬ 
tios. A Spearman rank correlation analysis is consistent 
with the null hypothesis that HCN/HCO’*' and Lir are 
uncorrelated {p = 0.021 ± 0.178, 2 ;-score= 0.057 ± 0.494; 
Table 3). The coefficients were computed using a Monte 
Carlo perturbation and bootstrapping method (with 10® 
iterations), as described by Curran (2014) and imple¬ 
mented in Curran (2015). 

The Spitzer/IRS observations of GOALS systems, pre¬ 
sented by Stierwalt et al. (2013), provide a calibrated 
measurement of the AGN contribution to the mid- 
infrared emission in each system, via the 6.2 p.m PAH 
EQW. In order to test claims that HCN (1-0) emis¬ 
sion is enhanced in systems with AGN, we explore the 
HCN/HCO"'' ratio as a function of this AGN diagnos¬ 
tic. In Figure 4 (Left), we plot HCN/HCO+ against the 
relative contribution of AGN and star formation to the 
infrared luminosity; where lower values of the PAH EQW 
indicate the presence of more energetically important 
AGN. We consider systems with PAH EQW < 0.2 pm to 
have energetically dominant AGN (Cmirmgn/^mir ^ 
60 %;) e.g., Petrie et al. 2011), while systems with PAH 
EQW > 0.55 pm are considered starburst-dominated 
(Brandi et al. 2006), and systems with intermediate ra¬ 
tios as composites, where Cmir is the mid-infrared lu¬ 
minosity and Ijmir,agn is the AGN contribution to the 
mid-infrared luminosity. 

We find that systems dominated by the AGN in the 
mid-infrared show elevated HCN/HCO'*' ratios, with a 
signal-to-noise ratio weighted mean ratio of 1.84 ± 0.43. 

and the HCN (1-0) and HCO”*" (1—0) (a 30% difference in beam 
size). 


PAH 6.2 fim EQW [/im] 
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Starburst dominated systems have a weighted mean ratio 
of 0.88 ± 0.28, while composite systems have a weighted 
mean ratio of 1.14±0.49. However, for systems which ap¬ 
pear to be star formation dominated, this ratio exhibits 
significant scatter. Several starburst and composite sys¬ 
tems have HCN/HCO“'" values which are comparable to 
the AGN dominated systems, suggesting that while en¬ 
ergetically dominant AGN are associated with elevated 
HCN (1“0) emission alone, the converse is not true: en¬ 
hanced HCN (1-0) emission does not imply the presence 
of an AGN. We note that although these starburst and 
composite sources with enhanced HCN (1-0) emission 
have substantial uncertainties in their HCN/HCO+ ra¬ 
tios, it is unlikely that the HCN/HCO+ ratio is simulta¬ 
neously overestimated for all six of these HCN-enhanced, 
composite/starburst sources. 

A Spearman rank correlation analysis shows 
HCN/HCO+ and the PAH EQW to be moderately 
anti-correlated, with p = —0.512 ±0.127 and a z-score of 
— 1.532 ±0.464 (Table 3). In Figure 4 (Right) we provide 
a “box plot” of the mean ratio, interquartile range, and 
full range of the HCN/HCO'^ values, separated by into 
source types based on the PAH EQW. The distributions 
of HCN/HCO’*' for pure starbursts and AGN-dominated 
systems show a clear offset, with the AGN-dominated 
systems show a relative enhancement of HCN (1-0) 
emission over the pure starbursts. 

It is worth noting there are fewer AGN-dominated 
sources than SB or composite sources; further observa¬ 
tions of low PAH EQW systems would be useful to ensure 
the existing objects are representative. 

4. WHAT DRIVES THE GLOBAL HCN (1-0)/HCO+ (1-0) 
RATIO? 

In Section 3 we demonstrated that globally enhanced 
HCN (1-0) emission (relative to HCO“'" (1-0)) is not 
correlated with the presence of an AGN (Figure 4). 
Costagliola et al. (2011) found similar results for a 
smaller sample of starbursts. Is there a straight-forward 
explanation for the observed global HCN (1-0) emis¬ 
sion in local (U)LlRGs? In the following subsections we 
explore a series of proposed explanations for enhanced 
HCN, including X-ray induced chemistry (Section 4.1), 
the presence of a compact, high-density source (Section 
4.2), and radiative pumping from absorption of mid- 
infrared photons (Section 4.3). We finish by discussing 
the possibility that the PAH EQW is not an ideal tracer 
of AGN (Section 4.4) and then mention future observa¬ 
tions which could be used to improve our understanding 
of HCN enhancements (Section 4.5). 

4.1. Comparison of HCN (1~0)/HCCC' (1-0) with 
X-ray Properties 

To consider the possible influence of XDRs resulting 
from powerful AGN, we investigated the HCN/HCO+ 
ratio as a function of X-ray properties, such as the hard¬ 
ness ratio and the total 0.5-10 keV X-ray luminosity, 
from Chandra X-ray observations of the GOALS sample 
(Iwasawa et al. 2011). The HCN/HCO+ ratio showed no 
correlation with either the hardness ratio or the X-ray 
luminosity, albeit with only ten sources common to both 
samples. Additional X-ray observations would be useful 
to further investigate the influence of XDRs. However, 
based on the currently available X-ray data, it does not 


appear that either the hardness of the X-ray spectrum or 
the total X-ray luminosity correlate with HCN/HCO’*', 
suggesting an XDR is not generally a major driver in en¬ 
hancing HCN (1-0), for these systems, possibly because 
the XDRs are spatially disconnected from the regions 
that dominate the global line luminosity. 

A complication is that enhancement from X-rays may 
be most effective in obscured AGN and diagnosing this 
activity is difficult with observations below 10 keV. Hard 
X-ray observations with NuSTAR (Harrison et al. 2013), 
though time-consuming, would provide an interesting 
test for the presence of obscured AGN in sources with 
enhanced HCN emission. 

4.2. HCN/HCO'^ Enhancements Through Source 
Compactness 

In compact environments HCO+ (1-0) appears to 
be more susceptible to self-absorption than HCN (1-0) 
(Aalto et al. 2015), which would lead to an increase in 
the observed HCN/HCO+ ratio. To test if elevated ratios 
are primarily associated with dense systems, we compare, 
in Figure 5, the HCN/HCO^ ratio with the [C ii]/Lfir 
ratios from the central Herschel spaxel (9".4 x 9".4) as 
measured by Diaz-Santos et al. (2013). In local galax¬ 
ies, the [C ii]/Lfir value decreases with increasing Lfir 
and dust temperature, Tdust- That the “[C ii] deficit” 
tracks Tdust suggests the deficit is the result of compact 
nuclear activity. Diaz-Santos et al. (2013) hnd this trend 
of a lower [C ii]/Lfir ratio with decreasing source size 
is present when considering only starbursts, suggesting 
that the [C ii] deficit is driven by the compactness of the 
starburst, rather than dust heated by an AGN. 

Therefore, if the HCN (1-0)/HCO''' (1-0) ratio is be¬ 
ing driven primarily by the density / compactness of the 
starburst, we should see a correlation between the two 
quantities. For sources with [C ii]/Lfir > 10“^, the 

signal-to-noise weighted mean to 

is 1.0 ± 0.4, while sources with [C ii]/Lfir < 10“^ 

(large [C ii] dehcits) have a mean ratio of 1.7 ± 0.5. For¬ 
mally, the ratio does not appear to vary as a function of 
[C ii]/Lfir though a Spearman rank analysis suggests a 
moderate anti-correlation (p = —0.401 ± 0.142, z-score= 
—1.165± 0.458; Table 3), though somewhat weaker than 
the anti-correlation of HCN/HCO’’' with PAH EQW. 
The majority of sources with [C ii]/Lfir < 10“^ have 
HCN/HCO’*' ratios > 1, consistent with a scenario in 
which a compact and dense starburst causes an enhance¬ 
ment of HCN (1-0). On the other hand, a substantial 
number of systems without significant [C ii]/Lfir deficits 
also have HCN/HCO'’' > 1. For composite and starburst 
systems (PAH EQW > 0.2 pm) higher HCN/HCO’*' ra¬ 
tios do not appear to be associated with lower [C ii]/Lfir 
values. Based on the available data, we cannot directly 
link the [C ii] deficit with enhanced HCN (1-0) emission. 

It is worth noting that HCN (1-0) can be enhanced 
in sources that do not show a strong [C ii] deficit; if 
a compact continuum source is surrounded by extended 
disk-like star formation, the system may have a “normal” 
[C ii]/Lfir ratio, but still show enhanced HCN (1-0) 
associated with the compact starburst on scales which 
cannot be resolved by Herschel. See Diaz-Santos et al. 
(2014) for a discussion on extended [C ii] emission in the 
GOALS sample. Assessing the spatial distribution of the 
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Figure 4. Left: {1-0 )/^hCO'^ (1 0) ^ function of the 6.2 PAH EQW from Stierwalt et al. (2013). The points are colored 

by their PAH EQW. The symbol shape notes the origin of the millimeter line measurements; circles are new data (see Table 2), diamonds 
are from Gracia-Carpio et al. (2006) and squares are from Costagliola et al. (2011). Right: Box plot showing the median (red line), 
interquartile range (boxes) and full range up to 1.5 x the interquartile range (IQR; black horizontal lines) for HCN/HCO"*” of the AGN 
dominated, composite, and starburst systems. A flier (a point with a value > 1.5 x IQR) is plotted with a ’+’ symbol. Upper and 
lower limits were not included in the box plot. HCN/HCO"*" is enhanced for systems which are AGN dominated in the mid-infrared, but 
some starburst dominated systems show similarly elevated ratios. Weighted by signal-to-noise ratio, the average ratio for AGN dominated 
systems (PAH EQW < 0.2 pm), pure starbursts (PAH EQW > 0.55 fim Brandi et al. 2006), and composite systems (0.2 fim < PAH 
EQW < 0.55 pm) is 1.84, 0.88, and 1.14, respectively. 



Figure 5. A comparison of HCN/HCO"^ with [C Ii]/Lpjp^, as 
a proxy for the compactness of the source. Systems with signif¬ 
icant [C ii] deflcits are located towards the left side of the plot. 
There is a dearth of sources with low [C Ii]/Lpjj:^ values and low 
HCN/HCO"^ ratios, but the scatter is significant and some systems 
with high [C Ii]/Lpjf^ values (less compact starbursts) still show 
elevated HCN/HCO"*" ratios, consistent with the mean of systems 
with substantial deficits ([C n]/Lpjp^ < 10“^). However there is 
significant scatter. Points are color-coded by the 6.2 pm PAH 
EQW. 

HCN (1-0)/HCO+ (1-0) would be a useful comparison 
with [C ii]/Lfir- do regions with “normal” [C ii]/Lfir 
correspond to regions with low HCN (1-0)/HCO“'' (1-0)? 

A possible concern with the above analysis is that the 
Lir, [C ii]/Lfir, and the PAH EQW are all correlated. 


Indeed, the Lir, the [C ii] deficit, and the 6.2 fim PAH 
EQW are interrelated, in that the most luminous sources 
tend to have the strongest [C ii] deficits (Diaz-Santos 
et al. 2013, because both Lir and [C ii]/Lfir correlate 
with the dust temperature;) and lower PAH EQWs (Pet¬ 
rie et al. 2011; Stierwalt et al. 2013). However, the con- 
I verse is not true; a normal [C ii]/Lfir value does not 

f guarantee a high PAH EQW. Systems with low or in- 

H termediate PAH EQWs are distributed across the range 

I of Lir and [C ii]/Lfir values seen for our sample (e.g., 

S Figure 5 and Diaz-Santos et al. 2013).. In contrast, 

§ essentially all the high PAH EQW (starburst) systems 

have normal [C ii]/Lfir and Lir ~ a few xlO^^ Lq. 
Thus, while these quantities are generally related, sys¬ 
tems with low or intermediate PAH equivalent widths 
appear to have a range of compactness (as diagnosed 
by [C ii]/Lfir) and so the PAH EQW and [C ii] deheit 
appear to be semi-independent. That the HCN (1- 
0)/HCO“*' (1-0) ratio does not correlate well with either, 
suggests the origin of HCN enhancement cannot be easily 
assigned solely to AGN or to the presence of a compact 
starburst. 

4.3. HCN/HCO'^ Enhancements Through Mid-infrared 

Pumping 

The strong mid-infrared continuum present in many 
(U)LIRGs may also influence these line ratios. The HCN 
molecule has degenerate bending modes in the infrared 
and can absorb mid-infrared photons (14 p,m) to its first 
vibrational state. The transitions have high level energies 
(> 1000 K) and mid-IR emission with a brightness tem¬ 
perature of at least 100 K is necessary to excite them 
(Aalto et al. 2007). A sign that infrared pumping of 
HCN is taking place is the presence of line emission from 
rotational transitions within the vibrational band (e.g., 
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Sakamoto et al. 2010). However, it may be possible for 
HCN to be infrared pumped without detectable emission 
from these rotational-vibrational lines. When the bright¬ 
ness temperature is close to the lower limit for pumping, 
the excitation of the rotational-vibrational line is so low 
that it takes a very large column density to result in a 
large enough optical depth for this rotational-vibrational 
line to be detected. 

Several sources in our sample do show evidence for in¬ 
frared absorption at 14 which could be associated 
with pumping of HCN (Lahuis et al. 2007; Sakamoto 
et al. 2010), but none of the starburst systems with en¬ 
hanced HCN emission show 14 fim absorption in the 
Spitzer/IRS observations from Inami et al. (2013). 

HCO"^ is similarly susceptible to infrared pumping, via 
a ~ 12 line; using the IRS high-resolution spec¬ 
troscopy from Inami et al. (2013) we searched for a 
^ 12 fj,m absorption feature in those systems with high 
ratios, where we postulated mid-infrared pumping could 
be important. We did not find any evidence for HCO'*' 
absorption. 

However, it is not clear that existing mid-infrared ob¬ 
servations are sensitive enough to establish firm upper 
limits on the importance of infrared pumping. Existing 
spectroscopy may not have the sensitivity to identify ab¬ 
sorption which is sufficient to affect level populations. 

4.4. Is the PAH EQW A Good Tracer of AGN Strength 

An obvious question is whether or not the PAH EQW 
is a robust measure of AGN strength in (U)LIRGs. Aalto 
et al. (2015) present interferometric observations of sev¬ 
eral objects, including the composite source CGCG 049- 
057, one of our sources with enhanced HGN emission. 
The HCN (3-2) V 2 = 0 and HCO+ (3-2) lines show ev¬ 
idence for significant self absorption and the HCN (3-2) 
V 2 = 1 ro-vibrational line is detected (consistent with 
mid-infrared pumping). Aalto et al. (2015) interpret this 
as evidence for a very compact (tens of pc) warm (> 100 
K) region surrounded by a large envelope of cooler, more 
diffuse gas. The total column density through these sys¬ 
tems is likely quite high (> 10^^ cm“^). If the column 
densities are sufficiently high to support mid-infrared 
pumping, the optical depths in the infrared may be high 
enough that the mid-infrared continuum level does not 
trace the intrinsic energy source - thus the PAH EQW 
diagnostic may miss highly embedded AGN. 

In Section 2.3.2 we note that the PAH EQW is an up¬ 
per limit to the AGN’s influence on large scales, but this 
measurement could also be considered as a lower limit to 
the AGN’s influence on smaller scales. Thus, consistent 
with results of Aalto et al. (2015) for GGCG 047-059, 
it is possible that these starburst or composite sources 
with high HCN (1-0)/HCO’'‘ (1-0) ratios host embed¬ 
ded AGN, and substantial mid-infrared optical depths 
result in an under-estimate of the AGN’s influence when 
using the PAH EQW. 

4.5. Future Observations 

Though we have shown in this study that 
some starburst-dominated systems have global 
Lrcn (i-o)/Lhco+ ( 1 - 0 ) excess of unity, 

the starburst-dominated systems outnumber AGN- 
dominated systems in the present sample. More 


IRAM 30m single-dish observations of AGN-dominated 
(U)LIRGs would be useful to study the dense-gas 
tracers in that population, particularly to see if 
the absence of AGN-dominated sources with low 
L(jcn (i-o)/^hco+ ( 1 - 0 ) truly representative or 

merely reflects small number statistics. 

Spatially resolved comparisons of the HGN emission 
with CO or HCO^ have been undertaken for systems 
known to host an AGN (e.g., Kohno et al. 2003; Iman- 
ishi et al. 2006, 2007, 2009; Davies et al. 2012; Imanishi & 
Nakanishi 2013, 2014); it would be instructive to perform 
the same exercise on starburst dominated systems. The 
Atacama Large Millimeter/sub-millimeter Array is the 
natural instrument for this, and such a study could inves¬ 
tigate the spatial variation in the 

ratio for starbursts. If the ratio peaks on the nucleus but 
is low elsewhere (as in systems with AGN), does that 
correspond to a compact mid-infrared emitting region? 
Spatially resolving this ratio in starburst systems with 
both high and low ratios may enable a separation of the 
relative importance of source compactness and/or mid- 
infrared pumping, in setting the ratio. 

Krips et al. (2008) studied the emission of multiple 
HCN and HCO+ transitions for twelve systems, find¬ 
ing evidence for systematic differences in the HCN (1- 
0)/HCO+ (1-0) and HCN (3-2)/HCO+ (3-2) ratios be¬ 
tween AGN and starbursts, with AGN having values > 2 
for both ratios. They also found evidence that the (3- 
2)/(l-0) ratio for both HCN and HCO®" varies between 
AGN and starbursts. The number of GOALS objects 
with (3-2) measurements of HCN and HCO®" is currently 
too small to draw conclusions regarding the (3-2)/(l-0) 
ratios as a function of AGN strength, but it would be 
possibly illuminating to compare the higher Jupper tran¬ 
sitions with the AGN diagnostic used here. 

The discovery of mid-infrared pumping and associated 
high-column densities in enhanced HCN sources sug¬ 
gests that even mid-infrared diagnostics such as the PAH 
EQW may miss highly embedded AGN. If this is the case, 
hard X-ray observations may provide the only conclusive 
evidence for AGN in sources with heavily obscured nu¬ 
clei. However, the weakness of the >10 keV X-rays in 
ULIRG AGN (e.g., Mrk 231; Teng et al. 2014) may make 
these observations difficult with existing facilities. 

Finally, the influence of infrared pumping is still uncer¬ 
tain. Future mid-infrared observations of these sources 
with The James Webb Space Telescope would provide 
tighter constraints on the importance and ubiquity of 
mid-infrared pumping, for both HCN and HCO®'. See 
also Aalto et al. (2007) for a discussion on using re¬ 
solved observations to distinguish mid-infrared pumping 
and XDR scenarios. 

5. THE RELATIONSHIPS BETWEEN MOLECULAR LINE 
LUMINOSITIES, STAR FORMATION RATES, AND GAS 
DEPLETION TIMESCALES 

5.1. Star Formation Rates 

We now turn to the relationship between the HCN (1- 
0) emission and the star formation rates in (U)LIRGs. 
While HCN (1-0) is generally taken to trace the mass of 
dense molecular gas, Mdense, previous studies have found 
conflicting results for the relationship between Lir and 
L(jcn ( 1 - 0 ) ■ Figure 6 (Left) we plot the standard re- 
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lation between L(jqj^ and Lir including additional 

measurements from Gao & Solomon (2004b, with val¬ 
ues corrected to our assumed cosmology and for the 3- 
attractor model of Mould et al. (2000)). 

We fit the relationship using a maximum likelihood 
technique, considering the errors in both Ljr and 
L(jcn ( 1 - 0 ) Hogg et al. 2010). When un¬ 

certainties for mm-line fluxes were not quoted in the lit¬ 
erature, we assume 20%. Uncertainties on the fit pa¬ 
rameters were determined using Markov Ghain Monte 
Carlo (MCMC) sampling^^ and quoted fit parameter un¬ 
certainties are for the 99% conhdence interval. 

To this point, we have not considered systematic un¬ 
certainties in our new observations, as the simultane¬ 
ous measurement of HCN and HCO’*" should result 
in signihcantly reduced systematic uncertainties in the 
HCN/HCO^ ratio,in the HCN/HCO“'" ratio, compared 
to non-simultaneous measurements. However, for the 
present comparison of Lir with these lines, the absolute 
flux of these lines may be subject to systematic uncer¬ 
tainties. We tested two approaches for dealing with sys¬ 
tematic uncertainties: 1) assigning additional systematic 
uncertainties (equal in magnitude to the statistical un¬ 
certainties) to each datapoint^^ and htting the relation¬ 
ship, or 2) including only statistical uncertainties^^ and 
allowing for additional uncertainties in the fitting pro¬ 
cess. In addition to fitting for the slope and intercept, 
we included an additional “nuisance” parameter, /, in 
the likelihood function to fit for the fractional amount 
by which the uncertainties are underestimated (e.g., due 
to not explicitly including systematic uncertainties). Af¬ 
ter the MGMG sampling, we marginalized over / when 
determining the final uncertainties for the slope and in¬ 
tercept. Both approaches yielded the same results for 
the best-fit relations, so we conclude that our fitting pro¬ 
cess appropriately accounts for non-statistical uncertain¬ 
ties and present the numerical results from the latter 
approach. For consistency, we show only the statistical 
uncertainties for values plotted in Figures 6 and 7. 

Our best fit relation to our new data combined with the 
literature data (omitting both AGN-dominated systems, 
in which Lir may not solely trace star formation, and 
HCN (1-0) upper limits) between Lir and is: 

lOgloLlR = (l-08lg j^g) lOgiQ -f (2.32^*^]^ Jg) 

( 1 ) 

This fit is shown in Figure 6 (Left) as the solid line. 
We also fit ^^_g^-LiR, hxing the slope to be lin¬ 

ear, which we show in Figure 6 (Left) as a dashed 
line. The slope of the (i_g)“LiR relation is con¬ 

sistent with being linear at the ~ 1.3(T level. We note 
that Garcfa-Burillo et al. (2012) found the (i-o)“ 

Lfir[40 — 400 fim] relation to be steeper than linear, 
with a slope of 1.23 ± 0.05 (68 % confidence interval). 

Calculated using the emcee python library (Foreman-Mackey 
et al. 2013). 

The values quoted by Garcia-Burillo et al. (2012) already 
include systematic uncertainties, so we did not increase the uncer¬ 
tainties on those points. 

In this case, we reduced the uncertainties on the Garcia- 
Burillo et al. (2012) points to remove their estimate for the sys¬ 
tematic uncertainties. 


In Figure 6 (Right) we perform the same compari¬ 
son, but utilizing instead of (i-o)' 

HCO’*' (1-0) has a higher critical density than CO (1- 
0), but lower than HCN (1-0), it is useful to test if it is 
also linearly tracked by the SFR. The (i-o)~^® 

relation is also consistent with linear: 

logio^iR = (0.94('(g j^g) logig Ljig;o+(i-o) + (3-5'^i2.2o) 

.( 2 ) 

Thus, our results are consistent with scenarios in which 
the HCN and HCO'*' emission both trace the dense gas 
associated with ongoing star formation, albeit with scat¬ 
ter, likely the result of one or more of the excitation 
mechanisms discussed in Section 4. 

5.2. The Lack of HCO^ Measurements for 
non-(U)LIRGs 

The inclusion of the Gao & Solomon (2004a) HCN ob¬ 
servations adds in a substantial number of galaxies with 
Lir < 10^^ L 0 , but corresponding HCO+ (1-0) obser¬ 
vations are not available. In order to assess the degree 
to which these lower-luminosity systems influence our fit¬ 
ting results, we also fit the HCN (I-O)-Lir relation with¬ 
out the points from Gao & Solomon (2004a); this ht is 
shown as the dot-dash line in Fig 7 (Left) and the best-fit 
relation is: 

logio Lir = (l-OSlg 3 ®) logig -|- (2.30('(3 04 ) 

( 3 ) 

which is consistent with that for the full dataset. 

5.3. Dense Gas Depletion Times 

The ratio L(jg,j 4 ^^ g^/LiR is often taken as oc 

Mdense/SFR, which corresponds to the depletion time 
ijdep) of the dense molecular gas. In Figure 7 (Left) 
we show L(jg,j 4 (- 4 _g)/LiR as a function of Lir, exclud¬ 
ing sources which are AGN dominated (PAH EQW 
< 0.2 fim). As in Figure 6 (Left), we include ob¬ 

servations from Gao & Solomon (2004a). The mean 
/ogi 0 (L(jQj^ (i_g)/LiR ) = I X 10“® with a RMS scatter 
of 0.22 dex. The AGN and starburst dominated sys¬ 
tems (as traced by the 6.2 fj,m PAH feature) have similar 
Lrcn ( 1 - 0 )/Lir ratios. 

Similarly, L(j^q+ ( 1 - 0 )/^^^ shows no obvious trend 
with Lir, but has substantial scatter. The mean 
logio(L'^^^+ ( 4 _g)/LiR ) = I X 10“® with a RMS scat¬ 
ter of 0.19 dex. 

The significant scatter in L(jg,]^ ^ 4 _Q^/L(jpQ+ can 

plausibly be attributed to the influence of multiple ex¬ 
citation mechanisms for these dense gas tracers, as de¬ 
scribed above. These observations are consistent with 
a scenario in which molecular abundance (Lepp & Dal- 
garno 1996), density (Meijerink et al. 2007), and radia¬ 
tive effects (Aalto et al. 1995) all influence the global 
HCN (1-0) emission, obscuring a simple link between 
Lrcn ( 1 - 0 ) ^Le mass of dense gas directly associated 
with ongoing star formation. Stated differently, the con¬ 
version factor between L(jg,j 4 ^^ g^ (or L(j^q+ ( 1 - 0 )^ 
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Figure 6. Ljj^ plotted as a function of (1-0) ^HCO'^ (1 0) points are color-coded by 6.2 pm EQW. The 

solid lines in the left and right panels denote the best fit relations of equations 1 and 2, respectively, while the dashed lines show the relation 
if we fix the slope to be linear. The fits do not include upper limits, but the (solid and dashed) fits to the (1-0) (left) 

includes the data from Gao & Solomon (2004b), corrected to our cosmology. In the left panel, we show a fit excluding the Gao Sz Solomon 
(2004b) points, as the dot-dash line. We find a slope consistent with being linear, when considering (^HCN (1-0)^ (Equation 1), 

consistent with Gao Sz Solomon (2004b). For Ljj^ ^^hCO"^ (1 0)^ ^ slope consistent with a linear relation (Equation 2). The fits 

do not include upper limits or AGN-dominated systems (PAH EQW < 0.2 pm). 
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Figure 7. Left: versus Ljj^. The Gao & Solomon (2004a) points have been corrected to the cosmology assumed here. 


Right: L' i , , /Ltr versus Ltr . Points are color-coded by the 

® HCO+ (1-0)' 
with Ljj^. 

Mdense likely depends on the relative HCN-H 2 abun¬ 
dance in addition to the overall density, excitation source 
(PDR vs XDR), and influence of infrared pumping. A 
better understanding of the processes in (U)LIRGs con¬ 
tributing to and are needed to 

determine if the scatter in Figure 7 is due to differences 
in the consumption rates of molecular gas or a vary¬ 
ing (^i_Q'^~Mdense convei'sion factor. Future obser¬ 

vations (Section 4.5), coupled with improved modeling, 
should be able to discriminate between these scenarios 


6.2 pm PAH EQW. We see no clear trend of the gas depletion time 

for the HCN and HCO+ emission in (U)LIRGs. 

6. SUMMARY 

We make use of new measurements of the putative 
high density gas tracers HGN (1-0) and HGO'*' (1-0) 
in a sample of local (U)LlRGs. A comparison between 
the ratios of these lines and the 6.2 fim PAH EQW 
mid-infrared AGN indicator suggests enhancements in 
the global HGN (1-0) emission (relative to HCO'*' (1- 
0)) does not uniquely trace the presence of an energeti¬ 
cally dominant AGN. While we find enhanced HGN (1-0) 























HCN (1-0) and HCO+ (1-0) in GOALS (U)LIRGs 


11 


emission relative to HCO+ (1-0) in objects hosting dom¬ 
inant AGN, we find the same magnitude of enhancement 
is also possible for systems which are dominated by star 
formation. The HGN (1-0) and HGO'*' (1-0) emission 
does not seem to be driven by a single process. It is likely 
that their emission is determined by the interplay of ra¬ 
diation field, gas column, and gas density. This hampers 
a simple interpretation of the line ratio. Existing data 
on the X-ray and mid-infrared properties of these sys¬ 
tems are not complete or deep enough (respectively) for 
us to prefer one of XDRs or mid-infrared pumping as the 
mechanism for enhancing the HGN (1-0) emission. 

We compare the HCN (1-0) emission with the star 
formation rate (Lm) and find a linear relationship, con¬ 
sistent with some previous studies (e.g., Solomon et al. 
1992; Gao & Solomon 2004b). However, our result is 
also consistent with a superlinear relationship, within our 
99% confidence interval, analogous to the recent study 
utilizing Lfir, by Garcfa-Burillo et al. (2012). The large 
scatter in the L'/Lir ratios is consistent with a scenario 
in which these dense gas tracers can be influenced by 
density effects, infrared pumping, and/or XDRs. This 
potentially complicates the determination of global dense 
gas masses and dense gas depletion times. 
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Table 1 

IRAM 30m Observing Log 


Source 

J2000 RA 
[deg] 

J2000 Dec 
[deg] 

z 

Obs Month 
JYYYY-MM] 

^int 

[min] 

T 

sys 

[K] 

Backend 

NGC 0034 

2.77729 

-12.10775 

0.0200 

2011-09 

276 

116 

WILMA 

MCG -02-01-052 

4.70904 

-10.36246 

0.0274 

2011-12 

127 

89 

FTS 

MGG -02-01-051 

4.71208 

-10.37672 

0.0270 

2011-12 

277 

90 

FTS 

IC 1623 

16.94804 

-17.50721 

0.0205 

2011-12 

95 

88 

FTS 

MGG -03-04-014 

17.53733 

-16.85272 

0.0342 

2011-12 

223 

86 

FTS 

IRAS 01364-1042 

24.72050 

-10.45317 

0.0487 

2011-09 

117 

98 

WILMA 

IC 214 

33.52279 

-1-5.17367 

0.0300 

2011-12 

244 

84 

FTS 

NGC 0958 

37.67858 

-2.93905 

0.0193 

2012-10 

61 

139 

FTS 

ESO 550-IG 025 

65.33333 

-18.81094 

0.0321 

2013-08 

298 

105 

WILMA 

UGC 03094 

68.89096 

-1-19.17172 

0.0247 

2012-10 

224 

126 

FTS 

NGC 1797 

76.93690 

-8.01911 

0.0149 

2012-10 

184 

113 

FTS 

VII Zw 031 

79.19333 

-1-79.67028 

0.0540 

2010-06 

42 

137 

WILMA 

IRAS F05189-2524 

80.25612 

-25.36261 

0.0435 

2011-12 

297 

107 

FTS 

IRAS F05187-1017 

80.27728 

-10.24619 

0.0283 

2011-12 

138 

88 

FTS 

IRAS F06076-2139 

92.44088 

-21.67325 

0.0376 

2013-08 

127 

107 

WILMA 

NGC 2341 

107.30000 

-1-20.60278 

0.0174 

2014-03 

61 

89 

FTS 

NGC 2342 

107.32525 

-1-20.63622 

0.0174 

2012-10 

326 

107 

FTS 

IRAS 07251-0248 

111.90646 

-2.91503 

0.0876 

2013-08 

255 

113 

WILMA 

NGC 2623 

129.60045 

-1-25.75461 

0.0185 

2010-06 

42 

200 

WILMA 

IRAS 09111-1007W 

138.40167 

-10.32500 

0.0564 

2013-08 

318 

103 

WILMA 

IRAS 09111-1007E 

138.41167 

-10.32231 

0.0550 

2011-09 

117 

98 

WILMA 

UGC 05101 

143.96500 

-1-61.35328 

0.0394 

2010-06 

191 

116 

WILMA 

CGCG 011-076 

170.30095 

-2.98396 

0.0247 

2012-10 

163 

123 

FTS 

IRAS F12224-0624 

186.26621 

-6.68103 

0.0257 

2012-10 

92 

135 

FTS 

CGCG 043-099 

195.46167 

-1-4.33333 

0.0374 

2011-09 

116 

98 

WILMA 

ESO 507-G 070 

195.71812 

-23.92158 

0.0215 

2011-12 

180 

110 

FTS 

NGC 5104 

200.34627 

-1-0.34248 

0.0186 

2012-10 

153 

121 

FTS 

IC 4280 

203.22212 

-24.20720 

0.0165 

2012-10 

132 

145 

FTS 

NGC 5257 

204.97042 

-1-0.84058 

0.0227 

2014-03 

132 

104 

FTS 

NGC 5258 

204.98854 

-1-0.82989 

0.0226 

2011-12 

180 

88 

FTS 

UGC 08739 

207.30800 

-1-35.25730 

0.0172 

2012-10 

102 

111 

FTS 

NGC 5331 

208.06750 

-1-2.10156 

0.0330 

2013-08 

286 

98 

WILMA 

CGCG 247-020 

214.93007 

-1-49.23666 

0.0258 

2014-03 

122 

153 

FTS 

IRAS F14348-1447 

219.40975 

-15.00633 

0.0830 

2011-09 

223 

139 

WILMA 

CGCG 049-057 

228.30455 

-1-7.22556 

0.0127 

2012-10 

71 

116 

FTS 

NGC 5936 

232.50360 

-1-12.98953 

0.0134 

2012-10 

91 

123 

FTS 

ARP 220 

233.73854 

-1-23.50323 

0.0181 

2010-06 

53 

159 

WILMA 

IRAS F16164-0746 

244.79913 

-7.90078 

0.0229 

2011-12 

201 

94 

FTS 

CGCG 052-037 

247.73545 

-1-4.08292 

0.0248 

2012-10 

214 

113 

FTS 

IRAS F16399-0937 

250.66754 

-9.72067 

0.0270 

2011-12 

191 

91 

FTS 

NGC 6285 

254.59998 

-1-58.95594 

0.0186 

2011-12 

127 

82 

FTS 

NGC 6286 

254.63146 

-1-58.93673 

0.0185 

2011-12 

53 

79 

FTS 

IRAS F17138-1017 

259.14900 

-10.34439 

0.0175 

2011-12 

149 

102 

FTS 

UGC 11041 

268.71599 

-1-34.77625 

0.0161 

2012-09 

41 

101 

FTS 

CGCG 141-034 

269.23598 

-1-24.01704 

0.0199 

2012-10 

61 

109 

FTS 

IRAS 18090-1-0130 

272.91004 

-1-1.52782 

0.0293 

2013-08 

286 

94 

WILMA 

NGC 6701 

280.80208 

-1-60.65312 

0.0132 

2012-10 

163 

165 

FTS 

NGC 6786 

287.72500 

-1-73.40992 

0.0250 

2011-09 

265 

95 

WILMA 

UGC 11415 

287.76833 

-1-73.42556 

0.0252 

2011-09 

287 

92 

WILMA 

ESO 593-IG 008 

288.62950 

-21.31897 

0.0487 

2011-09 

53 

120 

WILMA 

NGC 6907 

306.27750 

-24.80893 

0.0106 

2012-09 

102 

171 

FTS 

IRAS 21101-1-5810 

317.87667 

-1-58.38422 

0.0398 

2013-08 

308 

89 

WILMA 

ESO 602-G 025 

337.85621 

-19.03454 

0.0247 

2012-10 

142 

125 

FTS 

UGC 12150 

340.30077 

-1-34.24918 

0.0216 

2012-10 

81 

100 

FTS 

IRAS F22491-1808 

342.95567 

-17.87357 

0.0760 

2011-12 

158 

111 

FTS 

CGCG 453-062 

346.23565 

-1-19.55198 

0.0248 

2012-10 

122 

96 

FTS 

NGC 7591 

349.56777 

-1-6.58579 

0.0165 

2012-10 

101 

102 

FTS 

IRAS F23365-I-3604 

354.75542 

-1-36.35250 

0.0645 

2011-09 

32 

87 

WILMA 


Note. — Col 1: Source name, Col 2: Right Ascension, Col 3: Declination, Col 4: redshift, Col 
5: Year and month of observation. Col 6: Total on-source time, Col 7: System temperature, Col 
8: Backend used for measurements (either the Fourier Transform Spectrometer (FTS; Klein et al. 
2012) or the Wideband Line Multiple Autocorrelator (WILMA)). 
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Table 2 

IRAM 30m Fluxes 


Source 

S Tmbdv (CCH) 
[K km ] 

S T^bdv (HCN (1-0)) 
[K km ] 

S (HCO+ (1-0)) 

[K km ] 

S T^bdv (HNC (1-0)) 
[K km ] 

NGC 0034 

< 0.35 

0.71 ±0.12 

1.14 ±0.15 

< 0.44 

MCG -02-01-052 

< 0.33 

0.57 ±0.13 

< 0.27 

< 0.26 

MGG -02-01-051 

< 0.23 

0.17 ±0.05 

0.30 ±0.07 

< 0.16 

IG 1623 

1.40 ±0.14 

1.51 ±0.15 

4.00 ±0.15 

0.64 ±0.13 

MGG -03-04-014 

0.48 ±0.12 

0.94 ±0.12 

0.93 ± 0.09 

0.38 ±0.12 

IRAS 01364-1042 

< 0.46 

< 0.71 

< 0.45 

< 0.54 

IG 214 

< 0.15 

0.36 ±0.10 

0.56 ± 0.12 

< 0.14 

NGC 0958 

< 0.42 

< 0.71 

< 0.58 

< 0.69 

ESO 550-1 G025 

< 0.35 

0.45 ±0.11 

0.36 ± 0.11 

< 0.39 

UGC 03094 

0.63 ±0.12 

< 0.42 

0.74 ±0.14 

< 0.41 

NGC 1797 

< 0.37 

0.72 ±0.17 

0.55 ± 0.12 

< 0.43 

VH Zw 031 

< 1.00 

< 0.98 

1.26 ± 0.40 

< 0.96 

IRAS F05189-2524 

< 0.24 

0.73 ±0.13 

0.55 ± 0.13 

0.39 ±0.11 

IRAS F05187-1017 

0.49 ±0.12 

0.75 ±0.11 

0.76 ± 0.14 

< 0.18 

IRAS F06076-2139 

< 0.60 

1.05 ±0.20 

< 0.59 

< 0.50 

NGC 2341 

0.65 ±0.21 

< 0.35 

< 0.43 

< 0.60 

NGC 2342 

0.55 ±0.11 

0.68 ±0.11 

0.75 ± 0.13 

< 0.22 

IRAS 07251-0248 

< 0.41 

0.35 ±0.11 

< 0.33 

< 0.40 

NGC 2623 

< 1.27 

2.62 ± 0.48 

2.40 ±0.54 

< 1.23 

IRAS 09111-1007W 

< 0.36 

0.60 ± 0.08 

0.83 ± 0.16 

< 0.25 

IRAS 09111-1007E 

0.66 ± o.ie^- 

< 0.36 

< 0.41 

< 0.35 

UGC 05101 

0.99 ±0.20 

2.17 ±0.20 

1.25 ± 0.17 

0.92 ±0.20 

CGCG 011-076 

< 0.51 

1.09 ±0.15 

1.20 ±0.17 

< 0.35 

IRAS F12224-0624 

< 0.69 

< 0.55 

< 0.67 

< 0.66 

CGCG 043-099 

< 0.55 

< 0.54 

< 0.53 

< 0.53 

ESO 507-G 070 

< 0.48 

1.22 ±0.19 

1.87 ±0.21 

0.49 ±0.13 

NGC 5104 

< 0.67 

1.23 ±0.22 

0.74 ± 0.20 

< 0.41 

IG 4280 

0.55 ±0.15 

0.95 ±0.21 

1.17 ±0.23 

0.57 ±0.18 

NGC 5257 

< 0.29 

0.39 ±0.10 

0.46 ± 0.09 

0.37 ±0.11 

NGC 5258 

< 0.20 

0.43 ±0.10 

0.53 ± 0.07 

0.30 ±0.08 

UGC 08739 

< 0.51 

1.30 ±0.17 

0.94 ±0.19 

0.88 ±0.19 

NGC 5331 

< 0.31 

0.67 ±0.10 

0.45 ± 0.10 

< 0.43 

CGCG 247-020 

< 0.38 

0.96 ±0.15 

0.57 ±0.12 

0.35 ±0.09 

IRAS F14348-1447 

< 0.46 

< 0.52 

< 0.37 

< 0.51 

CGCG 049-057 

1.57 ±0.36 

3.21 ±0.32 

1.48 ± 0.20 

1.70 ±0.28 

NGC 5936 

< 0.48 

1.21 ±0.16 

0.81 ± 0.11 

< 0.46 

ARP 220 

4.37 ±0.62 

12.15 ±0.75 

6.02 ± 0.74 

7.85 ±0.60 

IRAS F16164-0746 

0.65 ±0.12 

0.54 ±0.11 

0.93 ± 0.12 

0.32 ±0.09 

CGCG 052-037 

0.57 ±0.12 

0.67 ±0.12 

0.74 ±0.10 

0.73 ±0.13 

IRAS F16399-0937 

0.66 ±0.14 

0.77 ±0.10 

0.68 ±0.10 

0.38 ±0.11 

NGC 6285 

0.38 ±0.10 

< 0.36 

< 0.25 

< 0.35 

NGC 6286 

0.93 ±0.20 

1.63 ±0.22 

1.67 ±0.21 

< 0.40 

IRAS F17138-1017 

1.01 ±0.17 

0.83 ±0.10 

1.03 ± 0.12 

0.47 ±0.11 

UGC 11041 

< 0.50 

1.51 ±0.26 

1.49 ± 0.20 

< 0.59 

CGCG 141-034 

< 0.51 

1.01 ±0.20 

< 0.49 

< 0.69 

IRAS 18090-1-0130 

< 0.30 

0.67 ±0.10 

0.78 ± 0.11 

0.28 ±0.08 

NGC 6701 

< 0.48 

2.77 ±0.22 

2.26 ± 0.22 

1.16 ±0.22 

NGC 6786 

0.30 ±0.10 

0.51 ±0.08 

0.56 ± 0.11 

0.35 ±0.09 

UGC 11415 

< 0.36 

< 0.28 

0.54 ±0.12 

< 0.22 

ESO 593-IG 008 

< 0.77 

< 0.98 

< 1.22 

< 1.05 

NGC 6907 

< 0.83 

2.29 ±0.31 

1.27 ±0.27 

< 0.65 

IRAS 21101-1-5810 

0.47 ±0.11 

0.53 ±0.09 

0.33 ± 0.08 

< 0.22 

ESO 602-G 025 

< 0.57 

0.84 ±0.21 

0.75 ± 0.12 

0.74 ±0.20 

UGC 12150 

0.69 ±0.15 

1.36 ±0.17 

1.59 ± 0.14 

0.51 ±0.14 

IRAS F22491-1808 

< 0.83 

< 0.82 

< 0.73 

< 0.71 

CGCG 453-062 

0.74 ±0.14 

0.97 ±0.16 

0.48 ±0.16 

< 0.42 

NGC 7591 

1.37 ±0.20 

1.35 ±0.18 

0.69 ± 0.13 

0.68 ±0.20 

IRAS F23365-I-3604 

< 0.58 

0.81 ±0.27 

1.28 ± 0.27 

< 0.88 


Note. — Col 1: Source name, Col 2: CCH (N = 1 —>■ 0) flux with Icr errors, Col 3: HCN (J = 1 ^ 0) flux with Icr 
errors. Col 4: HCO"*" (J = 1 ^ 0) flux with Icr errors. Col 5: HNC (J = 1 —)• 0) flux with Icr errors. The quoted upper 
limits are 3cr. 


^ Though formally a detection (> 3cr), this CCH flux appears to be spurious as the signal at the location of CCH is 
significantly broader than would be expected. 
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Table 3 

Spearman Rank Coefficients 


Quantities 

P 

2 -score 

^HCN (1-0)/^HCO+ (1-0) 
'"HCN (1-0)/'"hCO+ (1-0) 
'"HCN (1-0)/'"hCO+ (1-0) 

vs logi(,(LiR/L 0 ) (Figure 3) 

vs 6.2 (im PAH EQW (Figure 4, Left) 

vs logio([C iiJ/Lfir) (Figure 5) 

0.021 ±0.178 

-0.512 ±0.127 

-0.401 ± 0.142 

0.057 ± 0.494 

-1.532 ±0.464 

-1.165 ± 0.458 


Note. — Spearman rank coefficients for several relationships discussed in the text along with 2 - 
scores. Quantities and their 68 % confidence intervals were computed according to the Monte Carlo 
perturbation plus bootstrapping method (le5 iterations) discussed by Curran (2014) and using the 


associated code (Curran 2015). Upper and lower limits for (1-0 )/^hCO^ 


(1-0) 


were not included 


when computing these coefficients. 
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Figure 8 . Figure 1 Continued 
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Figure 9. Figure 1 Continued 
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Figure 10. Figure 1 Continued 
















































